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ABSTRACT 


The  algorithm  presented  Is  an  extension  of  the  Land  and 
Doig  branch  and  bound  method  combined  with  the  branch 
selection  techniques  presented  by  Beale  and  Small  to  solve 
integer  or  mixed-integer  linear  programs.  The  algorithm 
obtains  the  solution  by  solving  a  linear  program  with  upper 
and/or  lower  bounds  on  selected  branch  variables.  By 
systematically  changing  these  bounds,  and  maintaining  only 
the  current  canonical  form,  the  solution  is  assured  using  a 
minimum  of  excess  computer  storage  above  that  required  to 
solve  the  linear  programming  problem.  Thus  the  problem  can 
be  solved  entirely  within  the  computer  core,  and  the  problem 
converges  to  the  solution  faster  than  most  other  general 
integer  linear  programming  algorithms. 
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1.   INTRODUCTION 

In  the  literature  today,  there  are  many  algorithms  to 
solve  Integer  or  ml xed- Integer  linear  programs,  but  most  of 
these  algorithms  are  for  specialized  types  of  problems,  such 
as  transportation  problems  or  machine  scheduling  problems. 
As  Bellmore  and  Nemhauser  £  2  J  say,  "Integer  linear 
programming  algorithms  are  notorious  for  converging  rapidly 
on  one  problem  and  then  performing   miserably  on  the  next." 

The  algorithm  presented  here  Is  Harold  Oreenberg's 
natural  extension  of  the  Land  and  Dolg  branch  and  bound 
method  |J>J  '  CO  combined  with  branch  selection  techniques 
presented  by  E.M.L.  Beale  and  R.E.  Small  £tT|  ,  to  solve 
Integer  Linear  Programs.  In  addition,  the  algorithm  has 
been  developed  to  solve  mi  xed- Integer  linear  programs  where 
some,  but  not  all  of  the  variables  are    Integer-restricted. 

The  problem  is  solved  by  primal  linear  programming 
using  upper  and  lower  bounds  on  an  intelligently  selected 
integer-restricted  variable  that  does  not  have  an  Integer 
value  In  the  current  solution.  These  bounds  are  system- 
atically changed  until  an  optimal  solution  Is  assured.  In 
addition  to  speeding  convergence  by  having  less  branches 
than  most  other  methods,  this  method  also  uses  a  minimum  of 
excess  computer  storage  above  that  required  to  solve  the 
linear  programming  problem.  Thus  the  need  for  offline 
storage  Is  eliminated  and  the  entire  problem  Is  solved  and 
maintained  In  the  core,  giving  a  relatively  small  cpu  time. 


Computational  results  and  comparisons   with  other  algorithms 
are  included  to  show  this  small  time. 


11* FORMULATION  OF  THE  PROBLEM 

The  problem  to  be  solved  is:   find  X  ■  (x^/.../xn)  that 

minimizes    CX* 

(1) 
when        AX'=B 

0  1  Xj  1  Uj  ,   j-1,. ,.,n 

x.  Integer    for  some  or  all  j*l,...,n. 

C  is  an  n-vector,  B  is  an  m-vector,  and  A  is  an  mxn  matrix, 

and  the  u.  are  positive  integers  (or  are  infinite). 

To  formulate  a  problem  Into  this  form/  one  need  not 
restrict  the  slack/surplus  variables  to  being  integer/  thus 
there  is  no  need  to  restrict  the  elements  of  the  A,  B,  or  C 
vectors  to  being  integer. 

The  solution  to  (1)  includes  the  case  where  X;  is 
restricted  to  being  0  or  1.  This  accomplished  by 
restricting  x.  to  being  integer  with  an  upper  bound  of  one 
(i.e.  the  corresponding  u.  equals  one). 


III.   THE  BRANCHING  PROCEDURE 


A.   THE  LAND  AND  DOIG  METHOD 

The  branching  procedure  employed  In  this-  algorithm  Is 
an  extension  of  the  branch  and  bound  method  of  Land  and  Dolg 
|9|;  their  method  will  be  explained  first.  Land  and  Dolg 
first  solve  (1)  as  a  linear  program  to  obtain  the  continuous 
solution,  (z   ,  xi,...,xn).    If  the  continuous   solution  Is 

all  Integer,  the  problem  Is  solved.   If  there  are  fractions, 

o 
this  Initial  value  of  the  objective   function,  z  ,  serves  as 

a   lower   bound  on   the   value  of  any   feasible  optimal 

Integer-solution.   One  variable,  x»,   Is  then  selected  whose 

value  Is  not   Integer.   In  the  final   Integer-solution,  x:  , 

must   satlslfy  either  x.  <.  [x  •   ],  or   Xt  >.  [x  •  ]  ♦  1,  where 

[x .  ]  Is  the  greatest  Integer   less  than  x<  .   The  procedure 

Initially   Involves   the   solution  of   two   problems:    both 

contain  the  Initial  linear  programs,   but  with  an  additional 

constraint;   In  one   case  the  constraint  Is  x-.  ■   [x :  ],  and 

In      the  other      case      Itls      Xt      =      Qxjl      +      1.      Once      these 

solutions   are  obtained,   and  supposing   that   both  are   not 

Integer-solutions,  the  one  with   the  minimum  objective  value 

Is  Investigated  first,   while  the  other  Is   stored  for  later 

Investigation.    Suppose  the  minimum  of   the  two  objective 

functions  Is  z   ,  with  solution: 
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(z  1  ,x  }  ,  . . .,  [x  ?J  , . . .  ,x  ).   Now  a  new  variable,  x.  ,  whose 
a        j        n  k 

value  Is  not  Integer,  Is  selected,  and  three  new  problems 
are  solved:  these  are  the  Initial  linear  program  (1)  subject 
to  the  additional  constraints  x.  ■  [x >  ],  xk  =  [x^  } ,  In  the 
first  case,  x.  ■  [x?  J,  xk  *  Cx£  J  *  1,  In  the  second 
case,  and  x,  =  [x°  J  -  1,  in   the  third   case,  all  assuming 

The  procedure  Is  again  repeated  until  an  integer- 
solution  is  obtained.  At  every  step  where  there  Is  a  new 
variable  chosen  for  restriction,  there  Is  an  additional 
problem  to  solve.  Thus  the  next  step  would  have  four 
problems  to   solve,  and   the  one   after  that,   five,  If   the 

Integer-solution  Is  not  obtained.    Once  an  integer-solution 

q 
Is  achieved,   its  objective  function   value,  z   ,   becomes  a 

new  upper  bound   on  the  optimal  solution.    With  this  bound, 

the  procedure  then  investigates  the   other  problems  that  are 

solved  at   this  step.   At  any   point  in  the   branching,  when 

the    current   upper    bound   is    exceeded,    or   when    an 

Integer-solution  is  obtained,  that  branch  is  terminated. 

This  complete   process  enumerates  the   possibilities  for 

the  solution   to  (1)   In  a  directed   tree.   The   rooted  node 

corresponds  to  z   and  It  has   directed  branches  to  the  nodes 

corresponding  to   the  solution   of  the   linear  program   (1), 

with   Xj  set   at  a   specific   value,  each   of  these   Integer 

values   representing  a   single   branch.    At  each   of   these 

nodes,   the  branching   Is  repeated.    The  minimum  objective 

value  at  each  node  determines   the  next  branch,  which  limits 

the  potentially  large  size  of  the  tree. 
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The  problems  that  are  stored  for  further  Investigation 
rapidly  create  an  Immense  storage  problem,  thus  the  Land  and 
Dolg  method  becomes  unwieldy  for  problems  with  more  than 
Just  a  few  variables. 

The  algorithm  presented  here  takes  the  Initial  thought 
of  the  Land  and  Dolg  method,  but  uses  an  Inequality 
constraint  for  the  selected  variable.  This  Inequality  may 
be  as  an  upper  or  a  lower  bound  for  the  variable.  By 
employing  a  bounded  variable  linear  programming  solution 
technique  similar  to  that  from  Dantzlg  [  3  ],  the  algorithm 
allows  the  problem  to  be  solved  with  a  minimum  of  excess 
storage  used  above  that  storage  required  to  solve  the  linear 
programming  problem. 

B.   THE  NEW  BRANCHING  PROCEDURE 

As  with  the  Land  and  Dolg  method,  the  procedure  obtains 
the  continuous  solution  to  the  problem.  If  this  solution  is 
not  an  Integer  solution,  the  procedure  then  selects  a 
non-Integer  valued  variable,  x  j  ,  to  branch  upon.  This 
creates  two  new  problems;  both  contain  the  original  linear 
program  (1),  but  one  has  the  added  constraint  of  x  ? <.  [  x?], 
while  the  other  has  the  added  constraint  of  x.  2  C*?]  *  I* 
Initially,  only  one  of  these  problems  Is  solved.  The 
solution  to  this  problem  Is  treated  In  the  same  manner,  and 
the  process  Is  continued  until  an  Integer  solution  Is 
obtained,  or  the  problem  becomes  Infeaslble.  One  thing  to 
notice   Is   at   every   node   there   are  only   two   possible 
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branches,  and  there  are  no  problems  saved.  This  limiting  of 
branches,  alone,  tends  to  reduce  the  number  of  nodes 
required  for  the  solution  of  the  problem.  In  addition,  once 
an  Integer  solution  Is  obtained,  It  serves  as  an  upper  bound 
on  the  final  solution,  and  no  branching  continues  beyond  a 
node  where  the  objective  value  excedes  the  current  upper 
bound  solution. 

The  strategy  used  Is  to  Investigate  a  single  branch  of 
the  tree  until  zq  >.  zz  ,  where  zz  Is  the  current  upper 
bound  on  the  optimal  solution.  At  this  point,  backtracking 
begins,  to  find  a  node  that  has  only  one  previous  branching. 
This  backtracking  Is  accomplished  simply  by  removing  all 
bounds  from  the  nodes  that  are  further  out  on  the  tree,  and 
solving  the  resulting  problem.  This  solution  Is  the 
re-creation  of  a  previous  solution  that  would  have  been 
stored  otherwise.  The  second  branch  Is  then  made  utilizing 
the  Initial  problem  at  that  node  which  Is  re-created  from 
the  latest  transformation  of  the  equations,  and  the 
constraint  which  was  Initially  Ignored  at  that  node,  and 
solving  this  linear  program.  At  times,  the  procedure  may 
Impose  further  bounds  on  a  previously  bounded  variable. 
This  presents  no  problem  since  the  Initial  bounds  will  be 
satlslf led. 

Thus,  this  procedure  stores  only  the  current  solution 
and  the  current  bounds.  The  optimal  solution  Is  the 
existing  upper  bound  solution  when  It  becomes  impossible  to 
branch  any  further.  This  Is  realized  when  backtracking  to 
the  rooted  node  a  second  time  during  the  solution  procedure. 
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IV.   SIMPLEX  PROCEDURES  USED 

The  program  utilized  the  two  phase  simplex  method  to 
solve  (1) ,  together  with  a  bounded  variable  technique 
similar  to  that  In  [33.  By  using  this  technique,  the 
problem  maintains  m  constraint  equations  throughout  the 
solution  procedure  and  has  a  solution  where  some  of  the 
variables  will  be  at  their  upper  bounds,  some  of  the 
variables  will  be  basic,  and  some  of  the  variables  will  be 
at  their  lower  bounds. 

At  each  node  In  the  procedure,   we  are    Interested  In  the 

solution   to   (1)   which  satlslfles  the  new  bounds  given  by 

0  <.  1 .   <.  x.  <  u.  <   u  ,  1=1, ...,n,  where  u    Is  the  problem 

J      J  ""  J  "   J  -J 

upper  bound,  which  may  be  Infinite,   u.  Is  the  Imposed  upper 

bound  generated   by  the  branching   procedure,  and  1.   is  the 

Imposed   lower   bound,   also   generated   by   the    branching 

procedure.   After  achieving  a   solution   for  the   x.   which 

satlslfles  the  constraints  and  the   Imposed  bounds,   we  can 

obtain  the  following  canonical  form  [3],  from  the  constraint 

equations  In  (1): 


i   *  ^   yiixi 
1   j-m+1    J  J 


=  d.,  lsl,...,m, 


(3) 


-z    ♦  t.     c.  x.     ■  -z  , 

j-m+1   J   J 

where   the  x.  ,  l»l,...,m,  are   the  basic  variables,  and  the 

i 

x  ,  J»m*l,...,n,  are  the  non-basic  variables.   When  there  Is 
j 
a   feasible  solution   to  this   problem  (3),   with  upper   and 


H 


lower  bounds,   we  can  use  a   variation  of  the   usual  simplex 

method  to  obtain  optimal ity.   This  variation  is  suggested  by 

utilizing  a   method  to  improve   a  feasible  solution   that  is 

not   optimal.   To   improve  this   feasible   solution,  we   can 

increase  a  non-basic  variable  x    at   its  lower   bound   when 

J 
C   <  0,  or   we  can  decrease  a  non-basic  variable   x   at  its 

J  J 

upper   bound  when   c.   >  0,   where  c.   is   the  reduced   cost 

J  J 

coefficient  for  the   j tn  variable  as  found   In  the  objective 

function   in   the   cononlcal   form.    These   conditions   for 

optimal ity  come  from  the  following  theorem  |5|: 

Theorem: 

any  values  of  x  ,  }sl, , . , ,n,    satisifying  the  constraint 

J 
equations  in  (3)  and   0  <.  1  .  <.  x »  <.  u .   <.  u.,  j=l,...,n, 

are  an   optimal  solution  for   minimum  z,   with  objective 

vaiue  z  ,  Iff  c   >.  0  for  variables  at  their  lower  bound, 
o_      j 

1  ,  and  c .  <.  0  for  variables  at  their  upper  bound,  u  . 
The  proof  is  obvious,  for  any   Increase  In  a  variable  at  Its 
lower  bound,   or  any   decrease  In   a  variable   at  Its   upper 
bound    an  only   Increase  z.    More  detailed   proofs  may   be 
found  In  |3 | ,  or  In  | 6 | . 
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V.   THE  VARIABLE  SELECTION  PROCEDURE 

While  any  non-Integer  variable  In  the  present  solution 
can  be  chosen  for  the  branching  the  size  of  the  directed 
tree  can  be  further  reduced  by  selecting  the  one  which  will 
produce  the  optimal  solution  qulckest/  thus  reducing  the 
number  of  branchings  and  nodes.  The  method  employed  In  this 
algorithm  was  developed  by  Beale  and  Small  [U] ,  using  a 
suggestion  of  Drlebeck. 

The  process  Involves  comparisons  of  the  marginal  cost  of 
Increasing  and  of  decreasing  each  of  the  non-Integer  basic 
variables  In  the  current  solution,  as  In  the  cost  ranging 
procedure  of  linear  programming,  where  the  marginal  cost  of 
changing  each  of  the  non-basic  variables  Is  considered,  and 
the  non-basic  variable  with  the  greatest  change  which  still 
allow  feasibility  Is  entered  Into  the  basis.  The  marginal 
cost, J)  .,  of  Increasing  x.  is  given  by 

J>    =  Mln     (c.  /ly.-l), 
Ul  j,y   <  0    J     IJ 

U 

and  the  marginal  cost,!),,.,  of  decreasing  x.  Is  given  by 

■^  D  i  i 

J)  _.  =  Mln     (c.  /  y. .  ). 
D*  j,yjj  >  0   J     U 

where   c.  ,  is   the   reduced   cost   coeflclent  for  the   j 

variable  as  found  In  the  objective  function  In  the  canonical 

form.   By  applying  the  strategy  of  taking  the  variable  that 

makes  the  worst  of  the  two  alternatives  as  bad  as  possible, 

we  have  the  alternative  that  Is  the  most  unlikely  to  contain 

the  final   solution,  and  we  branch   on  this  variable   In  the 
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other  direction. 

This  strategy  Is  Impl  Imented  by  computing^    -p   ,xf. 

and  JL  ...    =J)irx  CI   -   f.),   for   each   non-Integer   basic 

variable  that  Is  restricted  to  being  Integer.    Where  f.  Is 

the  fractional  part  of  x.  ,  x.  =  Tx.  ]  *  f .  .   The  branching  Is 

ii   u  i J     i 

then  made  on  the  variable  that  makes  the  larger  of  these  two 

quantities  as  large  as  possible.   Say  that   Max^    > 

Max^   ,  we  then  go  back  to   find  I  such  that,L  =  Max  J^  #. 

Thus,  the  Ith  basic  variable  would  be  the  one  branched  upon, 

and  the  direction  of   the  branching  Is  x.  >  [x.    +  1,  since 

i      i 

It  Is  thej ,.  and  not  theJ     that  Is  the  greater. 

**"  Di  *■  Ui 
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VK   THE  ALGORITHM 

A.   THE  BRANCHING  PROCEDURE 

The  branching   procedure  Is   contained   In  the  following 
algorithm; 

1.  Definitions 

S(zz,x,  , . ,  ,,x  )  Is  the   current  upper  bound  solution 
1      n 

which  represents  a  feasible  solution  to  (1)  with  objective 
value  zz,  and  having  all  variables  at  Integer  values  which 
are  required  to  be  Integer  In  the  solution. 

I(!_)  Is  the  Index  of  the  L    bounded  variable. 

t  h 

B(L)  Is  the  bound  restricting  the  L  bounded 
variable. 

BSIGN(L)   Indicates   whether  the   bound   on  the   L 
bounded  variable  Is  an  upper  or  lower  bound:   BSIGN(L)  =  +  1 
for  an  upper  bound/  and  BSIGN(L)  =  -  1,  for  a  lower  bound. 

N(L)  Is  the  counter  to  show  the  number  of  branchings 
for  the  Lthbounded  variable;   Initially,  N(L)  =  0. 

2.  Logic  Procedure 

a.  Set  L=0,  and  take  zz  as  Infinite,  unless  there 
Is  some  known  feasible  Integer-solution  to  (1),  to  get  the 
Initial  S(zz,x-  ,...,x  ).  Solve  (1)  as  a  linear  program  to 
obtain  the  continuous  solution.  If  this  solution  has 
Integer  values  for  all  of  those  variables  that  are  Integer 
restr!cted/   the  problem   Is   solved.    Otherwise  go   to   b, 


18 


maintaining  the  canonical  form  of  the  solution  to  (1). 

b.  Set  L»L  ♦  1.   Go  to  b(l). 

(1)  Select  the   basic  variable,   x?  ,  that   Is 

J 

fractional,  but   required  to  be   Integer,  which  Is   the  most 
unlikely  to  produce  the  optimal  solution,  and  go  to  b(2). 

(2)  Set  BSIGN(L)  opposite  the  way  which  will 
produce  the  most  unlikely  chance  for  the  optimal  solution 
and  go  to  b(3)  or  b(U). 

(3)  If  BSIGN(L)  =  +  1,  set   B(L)  =   x?  .   Go  to 

J 
c. 

(«♦)   If  BS I  GNU)  «  -  1,  Set   B(L)  =   x?  ♦  1 .   Go 

J 

to  c. 

c.  Solve  the  linear  program  using  the  maintained 
canonical  form  with  the  new  bound  on  the  variable  x. , 
where  j  =  l(L).  Go  to  the  following  cases,  one  of  which 
will  hold: 

(1)  If  the  solution  produces  an  objective 
value,  z,  with  z  >  zz,  go  to  d. 

(2)  If  the  solution  produces  an  objective  value 

z,  with  z  <  zz,  and  the  solution  has  Integers  for  all  of  the 

required   variables,   redefine   S(zz,x„ , . . . ,x   )   as   a   new 

1      n 

feasible   Integer   solution   upper   bound,   where  zz  *  z,  x 

1 

,...,x  ,  Is  the  new  solution.   Go  to  d. 
n 

(3)  If  the  solution  produces  an  objective 
value  z,  with  z  <  zz,  and  the  solution  does  not  have 
Integers  for  all  of  the  Integer-restricted  variables,  go  to 
b. 
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(k)  If  the  problem  has  an  infeaslble  solution, 
go  to  d. 

d.   Set  LL  =  L,  and  to  to  d(l). 

(1)  If   N(L)=2,  go  to  d(2).    Otherwise  N(L)=1; 
go  to  d(3). 

(2)  If   L  =  l,   the   current   feasible   solution 

S(zz,x   ,...,x    )is   optimal,   and   the    algorithm   stops. 
1       n 

Otherwise,  set  L  =  L  -  1  and  go  to  d(l). 

(3)  Solve    the   linear    programming   problem 
starting  from  the   current  canonical   form   with   the  bounds 
B(l_),  B(L  +  1),...,B(LL),  removed.   If  j  =  l(L),  and  x.  does 
not  become  a  basic  variable,  the  solution  Is  non-unique;  x. 
Is  then  made  a  basic  variable.   Go  to  d(U)  or  d(5). 

(U)    If  BSIGN(L)   =   +1,  set   B(L)   =  B(L)   +1, 
BSIGN(L)  =  -1,  and  N(L)  =  2,  and  go  to  c. 

(5)    If  BSIGN(L)   =   -1,  set   B(L)   =  B(L)   -1, 
BSIGN(L)  =  +1,  and  N(L)  =  2,  then  go  to  c. 

This  completes  the  branching  part  of  the  algorithm. 

B.   THE  VARIABLE  SELECTION  PROCEDURE 

The  selection  procedure  for  which  basic  variable  to 
branch  upon,  and  which  way  to  branch  Is  contained  In  the 
following  algorithm. 

1.   Definitions 

I  Is  the  I    basic  variable. 

J  Is  the  j    column  In  the  basis. 

t  h 

IB(I)  Is  which  problem  variable  Is   the  I     basic 

variable. 
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IMU  Is  the  Index  for  which  I  has  the  max«£,  m  • 

IML  Is  the  Index  for  which  I  has  the  m.ax^_    . 

DU(J)  Js  the  current  ratio  ■  |c./y..|,lf  y, , <  0. 

DL(J)  Is  the  current  ratio  =  lc./yj.|,if  y{ . >  0. 

XDU  Is  the  current  value  of  the  mln  DU(J). 

J 

XDL  Is  the  current  value  of  the  mln  DL(J). 

XLDU  Is  the  current  value  of  the  maxj^ 

i     Ul 

XLDL  Is  the  current  value  of  the  max  1    _.  . 

i  ***  D  i 

BSIRN  Is  the  Indicator  for  whether  the  branch 
requires  an  upper  bound  or  a  lower  bound  on  the  variable  to 
be  restricted.  BSIGN  =  +1,  for  an  upper  bound,  and  BSIGN  ■ 
-1,  for  a  lower  bound. 

M  Is  the  number  of  constraint  equations/  I.e.  the 
size  of  the  basis. 


2.   Logic  Procedure 

a.  Set  XLDU=XLDL=IMU=IML=  0,  1=  0,  and  go  to  b. 

b.  Set  1=1  +  1.  If  IB(I)  is  an  artificial 
variable,  or  If  xiB/i\  's  not  restricted  to  being  Integer, 
or  If  x  ts  Integer,  go  back  to  a.  Otherwise,  set  J=0, 
XDU=XDL=  Infinity,  and  go  to  c. 

c.  Set  J  =  J  +  1.   If  c.  =  0,  repeat   step   c.   If 

J 

y..<  0,  go  to  c(l),  If  y..  >  0,  go  to  c(2),  otherwise,  y..  * 
0,  and  repeat  step  c. 

(1)  DU(J)  =  |c,/y..|.  If  DU(J)  <  XDU,  set  XDU 
■  DU(J),  and  go  to  c(3).  Otherwise,  go  straight  to  c(3), 
without  setting  these  values. 

(2)  DL(J)  »  ICj/yjjI.    If  DL(J)  <  XDL,  set  XDL 
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a  DL(J),   and  go  to  c(3).    Otherwise,  go  straight   to  c(3), 
without  setting  these  values. 

(3)  If  J  <  M,  go  to  c,  otherwise  this  Is  the 
last  comparison  for  the  I tn  basic  variable  and  go  to  d. 

d.  Go  to  the  following  cases. 

(1)  If  XDL  Is  still   equal  to  Infinity,  XLDL  Is 

unchanged,  otherwise,  If  XDL  x  f,R#n  >  XLDL,  set  XLDL  = 

XDL  x  f      and  set  IML  =  I,  go  to  d(2). 
I  B  (  I  ) 

(2)  If  XDU  is  still   equal  to  Infinity,  XLDU  Is 

unchanged,  otherwise,  If  XDU  x  (1  -  f|R(n>  >  XLDU,  set  XLDU 

=  XDU  x  (1  -fID,.  J/  and  set  IMU  =  I,  go  to  d(3). 
IBM) 

(3)  If  I  >  M,  go  to  b,  otherwise  this  Is  the 
last  basic  variable  to  consider,  and  go  to  e. 

e.  One  of  the  following  cases  will  hold. 

(1)  If  IMU  and  IML  both  equal  zero,  this 
procedure  was  called  In  error,  and  there  are  no  variables  to 
branch  upon. 

(2)  If  XLDU  >  XLDL,  then  the  variable  to  branch 

upon  Is  IMU,  and  BSIGN  =  +  1,  to  set  an  upper  bound  on  x    . 

I  MU 

(3)  If  XLDU  <  XLDL,  then  the  variable  to  branch 
upon  Is  IML,  and  BSIGN  =  -  1,  to  set  an  lower  bound  on  x,..,  . 

I  ML 

This   completes  the   variable  selection   portion  of   the 
al gor I thm. 
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VII.   EXAMPLE 

The  basic  branching  method,  along  with  the  selection 
procedure  for  which  variable  to  branch  upon  are  Illustrated 
In  the  following  ml xed- Integer  linear  programming  problem. 
This  Is  the  same  problem  used  In  [l]  to  Illustrate  the  Land 
and  Dolg  procedure. 

Mln    hx1    +  5x« 
when     3x.  +   x«  -  x,  =  2 

Xl  *  **x2      "  XU       =  5  *** 

3x.  +  2x?  -x5   ■  7 

x.  >.  0,  and  Integer,  I  =  1,...,5. 

The  continuous  solution  produces  the  canonical  form: 


x   +  -x   -  -x    =1& 
1   101*   10X5    10 


*2  "  lift  +  foX5   =10 


X3  +10X|*   10  5   "lT  (5) 


2   ♦  2   ♦U(    ..112 
lift    10  5      10 


Since  the  solution  has  fractions  for  at  least  one  of  the 
Integer  restricted  variables,  we  use  the  selection  procedure 
to  tell  which  variable  to  branch  upon,  and  which  way  to 
branch. 
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The  computations   to  select  the  branching   variable,  and 

the   direction   to  branch  are        Included   for   this   first 

branchl ng: 

1=1  f  »  0.8,   1  -  f  =  0.2. 

J  =  1 

c~  =  0,  thus  no  calculation. 
J 
J  =  2 

DL(2)  =  -g^   =  3.5, 

XDL    =  3.5. 

J  =  3 

DU(3)  -  hi     -2.75, 

XDU    =  2.75. 

XLDL  =  3.5  x  0.8  =  2.8;     IML  =  1. 

XLDU  =  2.75  x  0.2  =  0.55;   IMU  =  1. 

1=2  f  =  0.8,   1  -  f  =  0.2. 

J  =  1 

c   =  0,  thus  no  calculations. 
J 
J  =  2 

DU(2)  =  ixi  =  2.33, 
0.3 

XDU    =  2.33. 

J  =  3 

DL(2)  =  i-ti   =  ii.o, 

0.1 
XDL    =  11.0. 

XDL  x  f  »  11.0  x  0.8  =  8.8  >  XLDL. 

XDU  x  (1-f)  =  2.33  x  0.2  =  0.U67  <  XLDU. 
XLDL  =  8.8;     IML  =  2. 
XLDU  =  0.55;      IMU  =  1. 
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3  f  «  0.2,   1  -  f  =  0.8. 

J  =  1 


c.  =  0,  thus  no  calculations. 


J  =  2 

DL(3)  =  —   =  2.33, 
0.3 

XDL    =  2.33. 

J  =  3 

DUO)  =  i-ii   =  1.0, 

XDU    =  1.0. 

XDL  x  f  =  2.33  x  0.2  =  0.U67  <  XLDL. 

XDU  x  (1-f)  =  1.0  x  0.8  =  0.8  >  XLDU. 

XLDL  =  8.8;   IML  =  2. 

XLDU  =  0.8;   IMU  =  3. 

XLDL  =  8.8  >  0.8  *  XLDU, 

therefore:   BSIGN  =  -  1, 

IM  =  IML  «  2. 

Thus  we  Impose  the  lower  bound  on   x  >   1,  and  solve  the 

2~ 

problem  (5)  with  this  added   constraint  to  get  the  canonical 
form  for  node  1: 


2      1        7 

X   ♦  -x   -  -x    =   — 

13  2    3  5       3 


1        8 

x,  ♦  2x.  -  -x_   =-  —- 

«♦.     2    3  5      3 


x_  -   x2  "   x5   "   5  (6) 


7      k  28 

2   *  3*2  *  3*5   ="  "I' 
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with  (z,x1,x2,xz,xk,x5)    -  (35/3,5/3,1,11,2/3,0). 

Since  some  of  the  integer-restricted  variables  are  stili 

fractional,  we  again  use  the   selection  procedure  to  get  the 

new  bound  x,  2  2,  and  solve  the  new  problem   taken  from  (6) 

with  the  added  constraint,  x,  22,   to  get  the  canonical  form 

for  node  2: 

x5  -  2x2  -  3x,   s  -  7 

x   -  i»x_  -   x,   «  -  5 
li     2     1 

x,  -   x   -  3x,   »  -  2  (7) 

3      2      1 

-z  ♦  5x2  ♦  ifXj  *  0, 
where  (z,x.,x2,x5,x. ,x_  )  =  (13,2,1,5,1,1),  which  is  integer 
for  all  of  the  integer  restricted  variables,  thus,  this 
solution  becomes  the  solution  vector,  S(zz,x., • . .,x  ).  Next 
the  bound  x  2  2,  is  removed,  and  the  simplex  procedure  is 
applied  to  backtrack  to  the  previous  node  to  obtain  (6) 
again.  Now  we  impose  the  bound  Xj  ±  1,  and  solve  (6)  with 
this  new  constraint  to  get  the  canonical  form  for  node  3: 

xu  ♦   5xx  -  2x5  =   9 


x   ♦  2-x   -  ix    .  -Z 

2  2   12  5     2 

x   -  4*   -  zx    »  4  <8> 

3  2  1    2  5     2 

-z   -  g*t   *  |x5   -1U, 

where  (z,x  ,x  ,x  ,x  ,x  )  «  ( 1U,  1,2, 3,«»,  0) .   Since  z  >  13,  we 

again  backtrack  to  a  node  where  there  is  only  one  branch,  in 
this  case  to  node   0,  removing  the  bounds  x  <  1,   and  x  >  1/ 
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and   solving  by  the  simplex   procedure   to  again   get   the 

canonical  form   (5).   Then  a   new  constraint   of  x  £     0,  Is 

added  to  (6)  and  the  resulting   problem  is  solved  to  get  the 

canonical  form  for  node  k: 

x   +10x„  -  3x.   ■    8 
5      2      k 

x  ♦  kx     -   x   =    5 
1      2      4 

x   +llx   -   x    »   13  (9) 

3      2      k 

-z   -llx   ♦  Ux    ■-  20, 
2      d 

where  (z,x   ,x  ,x   ,x  ,x  )   ■  (20,5,0,13,0,8),   which   again 
1    2   3    i»   5 

has  z  >  13.    Since  there  are  no  nodes  with   just  one  branch 

which  can  be   backtracked  to,  the  problem  is   solved  and  the 

solution  is:    (zz,x  ,x   ,x  ,x   ,x  )   ■  (13,2,1,5,1,1).   The 

12345 

procedure  is  exhibited  in  the  tree  shown  in  figure  1. 

This  procedure  needed  only  two  nodes  to  obtain  the 
solution,  and  four  nodes  to  assure  the  solution,  while  this 
problem  (4),  as  shown  in  [2]  requires  five  nodes.  Thus  the 
advantages  of  this  procedure  over  the  Land  and  Doig 
procedure  is  evident  even  in  this  small  problem. 
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Node  4 

z  =  20 
z  >  zz 


Node  0 
z  »  11.2 


Node  0: 
Node  1 : 
Node  2: 
Node  3 : 
Node  4: 


Node  3 


z  =  14 
z  >  zz 


x2  >  1 


Node  1 
z  =  11.67 


Node  2 


z  =  13 
( Integer  Sol ) 
zz  =  z  =  13 


(112/10,18/10,8/10,42/10,0,0) 

(35/3,5/3,1,4,2/3,0) 

(13,2,1,5,1,1) 

(14,1,2,3,4,0) 

(20,5,0,13,0,8) 


Figure  1 


Had  we  just  used  some  arbitrary  rule  for  selecting  which 
non-Integer  basic  variable  to  branch  upon,  the  number  of 
nodes  may  have  been  much  greater.  The  following  tree,  In 
Figure  2,  Illustrates  the  arbitrary  rule  of  selecting  the 
first  non-Integer  variable  In  the  basis,  x  ,  and  then 
branching  above   this  variable   by   setting   the   bound, 


x   >  £x  1    +  1 . 

IT 
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Node  9 

z  =  Ik 
z  >  zz 


Node  0 
z  -  11.2 


Node  2 
z  =  11.75 


Node  8 

I nfeasible 


Node  7 

z  =  13 
( I nteger  Sol ) 
zz  ■  z  »  13 


Node  6 


z  =  20 
z  >  zz 


Node  5 


z  =  17 
( I nteger  Sol ) 
zz  a  z  =  17 


Figure  2 
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VIM.   COMPUTATIONAL  RESULTS 

The  algorithm  has  been  programmed  In  Fortran  I V,  and  run 
on  an  IBM  360/67.  The  complete  program  is  given  in  the 
appendix.  The  computer  core  needed  to  run  a  problem  is  62K 
plus  36n  ♦  32m  ♦  km  x  (m  +  n)  bytes,  where  n  is  the  number 
of  variables,  and  m  is  the  number  of  constraints  in  the 
particular  problem  when  in  the  form  as  (1).  The  program 
utilizes  a  simplex  procedure  in  real  variables  with  a 
round-off  correction  which  requires   that  all  problems  using 

this  program  have  all  integer  coefficients  when  in  the  form 
as  (1). 

The  results  in  Figure  3  are  for  problems  in  [7],  and  the 
data  for  the  other  programs  are  from   [7]  and  f.8] . 

The  results  show  that  this  branch  and  bound  algorithm  is 
faster  than  most  of  the  general  algorithms  for  integer  or 
mi xed- integer  linear  programs. 
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Comparison  of  Computational  Experience  on  Haldi  and  IBM 


Source 

m 

n 

LIP1 

IPM2 

IPM3 

BBA 

BVA 

(g> 

(h) 

(a) 

(b) 

(c) 

(d) 

(e) 

1  ter 

Iter 

1  ter 

1  ter 

Time 

Time( j ) 

(f) 

(sec) 

(sec) 

Haldi 

1 

k 

5 

25 

19 

51 

91+ 

3.15 

2 

k 

5 

21 

18 

58 

26 

.83 

3 

k 

5 

22 

17 

80 

kl 

1.33 

— 

k 

k 

5 

17 

29 

35 

12 

.38 

— 

5 

6 

5 

181 

765 

F(i) 

22 

30 

6 

6 

5 

126 

20 

F 

2U5 

8.2 

12 

7 

k 

5 

173 

731 

F 

31 

1.03 

21* 

8 

k 

5 

125 

20 

F 

31 

1.02 

9 

9 

6 

6 

U3 

39 

192 

22 

8 

10 

10 

12 

120 

F 

F 

17 

11 

21 

56 

101* 

-- 

— 

-- 

12 

21 

56 

103 

— 

— 

— 

13 

21 

56 

683 

— 

— 

— 

lit 

21 

56 

106 

-- 

— 

— 

15 

21 

56 

877 

— 

— 

-- 

IBM 

1 

7 

7 

11 

7 

8 

18 

2 

7 

7 

30 

10 

17 

21 

3 

3 

l» 

63 

2 

13 

1*8 

1.7 

k 

k 

15 

15 

60 

2k 

2k 

23 

5 

15 

15 

332 

l»53 

1078 

151* 

6 

31 

31 

3«*1 

60 

1*17 

-- 

7 

12 

50 

391 

F 

7k 

-- 

8 

12 

37 

63 

1*53 

33 

-- 

9 

50 

15 

957 

5561 

F 

— — 

in  this  paper, 
in  [8j,  times  are  from  T8]  , 
pivot  steps  required  to 


(a)  Values  are  from  [7] . 

(b)  Values  are  from  [71. 

(c)  Values  are  from  L7J  . 

(d)  The  algorithm  presented 

(e)  The  algorithm  discussed 

(f)  Iterations  are  the  total 
solve  the  problem. 

(g)  The  number  of  constraints. 

(h)  The  number  of  variables,  excluding  slacks.  Each  prob- 
lem has  as  many  slacks  as  it  does  constraint  equations 

(i)  An  "FM  means  that  the  problem  was  unsolved  after  3000* 
I terat ions. 

(j)   Times  to  solve  problem  on  IBM  360  series  computers. 


Figure  3 


31 


CCMMON  A(  25,  25), B(  25,  25),IB(  25), H(  25),AX(  25) 
COMMON  X(  25),XLB(  25),XUB(  25),XI(  25), Y(  25), R(  25) 
COMMON   XSOL(  25),IRV(  25) 

COMMON  MM,N,M, JS , E , I R , AD, W, Z, IM,DI, ITR ,XH, AH, HH, EE , BND 
COMMON  ZZ, DJS, ITER.L, I W,NODEO, I NTRE S, XT, XK , I D ,RL , BS IGN 
CCMMON   IREAD, IWRIT E, IBNB 
DIMENSION  ITITLE(15) 
DATA  ISTOP/'STOP1/ 

C  PROGRAM  TO  ENTER  DATA  AND  CALL  LP. 

C  MM  =  NUMBER  CF  CONSTRAINT  EQUATIONS. 

C  N=NUMBER  OF  VARIABLES — THE  #  OF  ELEMENTS  IN  THE  X-VECTOR. 

C  BND  =  0  MEANS  NO  BOUNDS. 

C  3ND  =  1  MEANS  BOUNDS. 

C  IL  IS  A  CONTROL  VARIABLE  TO  SIGNAL  THE  END  OF  DATA  FOR  A 

C  PARTICLUAR  READ  STATEMENT.   IL=1  FOR  LAST  CARD,  0  OTHER. 

C  INTRES  =  C  MEANS  THAT  THE  PROBLEM  IS  A  ALL  INTEGER 

C  RESTRICTED  PROBLEM. 

C  INTRES  =  1  MEANS  THAT  THE  PROBLEM  IS  A  MI XED- I NTEGER 

C  PROBLEM. 

C  IRV(I)  =  0  MEANS  THAT  THE  I-TH  VARIABLE  IS  NOT  INTEGER 

C  RE  STR I CTED • 

C  IRV(I)  =  VMEANS  THAT  THE  I-TH  VARIABLE  IS  INTEGER 

C  RESTRICTED. 

C  NOTE  THAT  THE  SLACK/SURPLUS  VARIABLES  DO  NOT  HAVE  TO  BE 

C  INTEGER  RESTRICTED. 

C  WHEN  INTRFS=  0,  THE  VALUE  OF  IRV(I)  IS  IMMATERIAL,  THE 

C  PROGRAM  WILL  HAVE  ALL  IRV(I)=0,  BUT  WILL  DISREGARD  THIS. 

C  M  IS  THE  VARIABLE  TO  AUGMENT  THE  MARTIX  FOR  THE  OBJ.  FUNC. 

C  H(I)  ARE  THE  RIGHT  HAND  SIDES, HIM)  =  INITIAL  Z. 

C  H(M)=INITIAL  Z,  WHICH  IS  ZERO  UNLESS  OTHERWISE  PUT  IN  PRO. 

C  A(M,J)  ARE  THE  COSTS. 

IWRITE  =  8 
IREAD   =  4 

2  FORMAT  (2I5,F5.0, I5,15A4) 

4  FORMAT( 11,213, F5.0) 

6  FORMAT  (  II, 13, F5. 0,12) 
10  FORMAT! 1H  ,10F8.0) 

14  F0RMAT(1H0,10X,  •  HAVE  FINISHED  •/•1,» 
49  FORMAT  ( 1H 1 , 3 5X, 1 5A4// ) 

51  FORMAT  ( •  MINIMIZE:* ) 

52  FORMAT  (•  SUBJECT  T0:») 

53  FORMAT  ( 1X,24F5.0) 

54  FORMAT  ( • + ■ , 124X , ■ = • , F5.0 ) 

55  FORMAT  (•  THIS  PROBLEM  HAS', 15,'  CONSTRAINT  EQUATIONS, 
1  ,15, •  VARIABLES,  AND  NO  BOUNDS*,/) 

56  FCRMAT  (•  THIS  PROBLEM  HAS', 15,'  CONSTRAINT  EQUATIONS, 
1  ,15, •  VARIABLES,  AND  HAS  UPPER  BOUNDS',/) 

57  FORMAT*  •  WITH  THE  FOLLOWING  UPPER  BOUNDS  ON  THE  VARP, 
l'BLES  (A  -1  INDICATES  NO  BOUND  ON  THAT  VARIABLE):1) 

58  FORMAT  I •  THE  FOLLOWING  VARIABLES  MARKED  WITH  A  ONE  A1, 
l'PE  INTEGER  RESTRICTED:') 

59  FORMAT  (IX, 2415) 

15  READ  (IREAD, 2)  MM , N, BND, I NTRE S , I T I TLE 
IF  (  ITITLE(l) .EQ.ISTOP)   STOP 

M  =  MM  +  1 

DO  1  1=1, M 

H(  I)=0. 

DO  1  J=1,N 
1  A( I,J)=C. 

DO  9  J=1,N 

IRV(J)  =  0 

XUB( J)=-l. 
9  XLB(J)=C. 

3  READ     (IREAD, 4)  IL,I,J,A(I, J) 
IF  { IL.  EQ.OGO    TO    3 

5  READ    (IREAC6)  IL  ,  I,H(  I  ) 
IF  ( I  L.  EQ.OGO    TO    5 
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7  READ    <IREAC,6)     IL  ,  J  ,  A  (  M,  J  )  ,  I  RV  (  J  ) 
IFIIL.EQ.DGO    TO    7 
IFCBND.EQ.OGO    TO    8 

11    READUREADf6  )     IL,J,XUB<J) 
IF ( IL.EO.C )GO    TO    11 

8  WRITE     (IWRITE,49)     ITITLE 

IF    (BND.EQ.O.)     WRITE    (IWRITE»55)    MM,N 

IP     IBND.EQ.l.)     WRITE     (IWRITE,56)     MM,N 

WRITE    <IWPITE,51) 

WRITE  (IWRITEf53)  ( A ( M, J ) , J= 1 , N  ) 

IF  (BNO.EO.O)  GO  TO  61 

WRITE(  IWRITE,57) 

WRITE     (IWRITE,53)        ( XUB (J ) , J=l , N ) 

61  IF    ( INTRES.EQ.O)       GO    TO    62 
WRITE    UWPITE.58) 

WRITE     (IWRITE,59)  (  IRV ( I) , 1  =  1 , N ) 

62  WRITE(  IWP  ITE, 52) 
DO  6C  L  =  liMM 

WRITE  <IWRITE,53)  < A< L  .  J  )  , J  =  l , N ) 
60  WRITE  (  IWR  ITE, 54)   H(L) 
C  EXECUTION  TIMED  FROM  THIS  POINT 

CALL  LP 

WRITE  <IWkITE,14) 
C  END  OF  EXECUTION  TIMING 

WRITE  (IWRITE,A9)  ITITLE 

DO  47  J=1,N 

IFtXSOLC J) )48, 47,48 
48    WRITE(  IWRITE,50) J,XSOL( J) 
50    FCRMAT(lHr,5X,5HXSCL( ,I2,2H)=,F13.4) 
47    CONTINUE 

WFITE  ( IWRITE,45)ZZ 
45  FCRMAT(1HC,3HZZ=,F10.5) 

GC  TO  15 

END 

SUBROUTINE  LP 

COMMON  A(  25,  25), B(  25,  25),IB(  25), H(  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB(  25),XI(  25), Y(  25), R{  25) 
COMMON   XSOL(  25),IRV(  25) 

COMMON  MM,N,M,JS,E,IR,AD,W,Z,IM,DI , I TR ,XH , AH ,HH, EE , BND 

COMMON  ZZ,DJS, ITER,L, IW,NODEO, INTRES , XT , XK , I D, RL , BSIGN 
COMMON   IPEAD, IWRITE,IBNB 

C  THIS  IS  THE  SUBPROGRAM  WHICH  CONTROLS  THE  SIMPLEX 
C  PROCEDURES,  THE  INITIAL  SET-UP  OF  THE  COMPUTATIONAL 
C  ARRAYS,  ANC  THE  TERMINATION  OF  THE  PROBLEM. 

NODEO  =  0 

ZZ=10C0000. 

IW=0 

Z=H(M) 

L=0 

M=MM+2 

ID=C 

DO  3  J=1,MM 

DO  2  1=1, M 

2  B(  I,J)=C. 
B(M, J)=-l. 

3  B(J,J)=1. 

C  XI(J)=*1.  MFANS  INCREASE. 
C  XI(J)=-1.  MEANS  DECREASE. 
C  XI (J)=2.   MEANS  BASIC. 
C  XI(J)=0.  MEANS  IMONBASIC. 
DO  24  J=1,N 
XI(J)=0. 
24  X(J)=0. 

MM 


DO  4  1=1. 
Y(I)=H(I  ) 
IB  < I )  =  -I 
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6 

C  W  = 

14 


8 
5100 


15 

INCR 
19 

16 

OECR 

17 
61 


10 
9 


213 
221 


3  6 
39 
40 
38 

62 
28 
44 
50 


AD  =  1. 

DI  =  1. 

E=.5 

W=0. 

DO  5  J 

A(M,J) 

DC  5  I 

A(M, J) 

DC  6  I 

W=h+H( 

SUM  OF 

WRITE( 

FORMAT 

ITR  =  C 

ITER=0 

JS=0 

IF  (  IT 

FORMAT 

IF  (  IT 

DO  9  J 

IF ( XI ( 

D=0. 

DC  7  K 

D=D+B< 

IF(W.E 

IF(ABS 

IF(D)1 

IF(XUB 

IF(ABS 

EASE. 

XI(J)= 

GO  TO 

IFiABS 

EASE. 

XI  (J)  = 

IF(JS. 

JS  =  J 

DJS=D* 

GO  TO 

IF1DJS 

CCNTIN 

IF(JS. 

WPITE( 

FORMAT 

CALL  E 

ITER -I 

IF  (L. 

IF(IR. 

IF (BND 

IF  (  BN 

IF(XH. 

DO  38 

IF(IB( 

Y(  I)=Y 

GO  TO 

K=IB(I 

X(K)=X 

CONTIN 

X( JS)= 

IF(W.E 

W=W+AX 

Z=Z+AX 

WRITE( 

FORMAT 

IF  I IR) 

CALL  P 

CALL  F 

IFUR. 

IF(W.G 

IFUW. 

w=o. 


=  1,N 

=0. 

=  1,MM 

=A(M,J)-A< I , J) 

=  1,MM 

I) 

THE  H(I  ).  1  =  1, 
IWRITE,14)W 
<1H0,5X,  •  INITIAL 


..,MM. 

W=  •  ,F12.4) 


ER.GT.225)  WRITE  ( I WRI TE f 5100) 

(1CX  ,•  ITERATIONS  ARE  ABOVE  LIMIT  SET  IN  PROG.1) 

EP.GT.225)  RETURN 

=  1,N 

J) .EQ.2.  I  GO  TO  9 

=  1,MM 

M,K)*A(K,J) 

CO.  )  D=D+A(M,J) 

(D) .LE.E)GO  TO  9 

5,9,16 

(J).E<i).-l.  )GO  TO  19 

(XUB(J)-X(J  )  ).LT.E)GO  TO  9 


1. 
17 
(X{J  ) 


XLB( J) ).LT.E)GO  TO  9 


-1. 
GT.C  )GO 


TO  10 


XI(  J 
q 

-D*X 
UE 

EQ.C 
I  WRI 
(1H 
LEM 
TEP  + 
LT.C 
GT.C 
.EQ. 
D.EQ 
EO.l 
C 
1=1, 
I)  )3 
(I)- 
38 
) 

(K)- 
UE 

X(  JS 
Q.C. 
(M)* 
(MM+ 
IWRI 
(1H 
528, 
IVOT 
LAG 
EQ.O 
T.E) 
EQ.l 


K J )  .GE.EJGO  TO  61 

)GO  TO  101 

TE,213) JS,DJS,XI( JS) 

,5X,3HJS=,I5,5X,4HDJS=,F1D.4,5X,3HXI  =  ,F1C4) 

1 

)  GO  TO  529 

)GO  TO  36 

CIGO  TO  33 

. 1. AND. L. EQ.O)   GO  TO  33 

000000. )GO  TO  33 

ALL  FLAG  (1105) 

MM 

9,528,40 

AX( I >*XH*XI< JS) 


AX(  I  )*XH*XI ( JS) 

)+XH*XI ( JS) 
)G0  TO  62 
XH*XI(JS) 

i)*xh*xh  js) 

te,28)w,z,iter 

,10x,2hw=,f10.4,10x,2hz=,f10.4,5x,5hiter=,i5 

50,44 

(  3083   ) 
)  XI(JS)=0.0 
GO  TO  8 
)G0  TO  8 
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M=M-1 
IW=1 

WRITF(  I  WRITE, 73) 
73  FORMAT(1HO,10X,  •  START  PHASE  TWO  »l 
GC  TO  8 

33  WRITE( IWRITE,704) 

WRITE  (6,5701) L,R(L),Z,ZZ 
704  FGRMAT(1HG,10X,18HUNB0UNDED  SOLUTION) 

RETURN 
101  IF(W.GT.E)GO  TO  35 

WRITE( IWRI TE,204) 

204  FORMAT ( 1HC ,5 X ,6HRESULT ) 
DC  31  J=1,N 

IF(ABS(X( J) ).LT.E)GO  TO  31 
WRITE!  IWRITE,200) J,X( J) 
200  FCRMATUH  ,  5X  ,2HX  (  ,  I  2  ,  2H  )  =  ,  F10  .4  ) 
31  CONTINUE 

Z=SlGMAINT(DI*ABS(Z)  +  .5)/DI,Z) 
WRITE( IWRITE,527)Z 
527  F0RMAT(1H0,5X,2HZ=,F17.5) 
;  HAVE  CONTINUOUS  SOLUTION. 

IF  (  INTRES.EQ.  1)  GO  TO  301 

IF  (Z+1-E.GE.ZZ)  GO  TO  134 

301  IF  (Z+E.GE. ZZ)    GO  TO  134 

IFCDI-1. 1601,601,603 
601  WRITE(  [WRITE, 7011 

WRITE  (6,5701) L,R(L),Z,ZZ 
5701  FORMAT  ( 5X  ,  • R(  • , I  3, • )  =  • , F 1 4. 6 , 5X , ■ Z  =  LF11.4, 
15X,*ZZ  =  LF11.4,/) 
7C1  FORMAT  ( 1HC, 5X , 22HHAVE  INTEGERS  REQUIRED) 

GC  TO  51 
603  IF  ( INTRES.EQ. 1)  GO    TO  550 

IF  (ABS( ABS(Z)-AINT(ABS(Z)+E) ) .GE.E)  GO  TO  6032 
550  DO  637  1  =  1, MM 

IF  (  INTRES.EQ.  LAND.  IRV(  IB(  I)  ).EQ.O)  GO  TO  637 
IF (  IB(  I) )637,637,83 
83  J=IB(I ) 

IF(ABS(X( J)-AINT( ABS(X(J) )+E) ) -E ) 63  7, 637, 6032 
637  CONTINUE 
;  HAVE  ALL  INTEGERS  REQUIRED. 
WRITE( IWRITE,701) 
WRITE  (6,5701)L,R(L) ,Z,ZZ 
51  DO  45  J=1,N 

XSOL(J)  =  AINT  (X(J)  +  0.5) 
45  IF  (  INTRES.EQ. LAND. IRV( J) .EQ.O)   XSOL(J)  =  X(J) 
ZZ  =  Z 

IF  (BND.EQ.O)   GO  TO  46 
IF  (BND.EQ.l. AND. L. EQ.O)  GO  TO  46 
GO  TO  34 
6032  WRITE(  IWPITE,703) 

WRITE    (6,5701)L,R(L),Z,ZZ 
703    FORMATS     HAVE    FRACTIONS1) 
BND=1. 

IFCID.EQ.500 )  GO  TO  41 
42  CALL  BOUND 

GO    TO    8 
35    WRITE(  IWRITE,205) 

WRITF  (6,5701) L,R(L),Z,ZZ 

205  FORMAT (1H0,»  INFEASIBLE*) 
IF(BNC.EQ.O. )G0  TO  528 

GO  TO  34 
134  WRITE  (  IWRITE.702) 

WRITE  (6,5701)  L,R(L),Z,ZZ 
702  FORMAT  (•  Z  VALUE  GREATER  THAN  SOLUTION  UPPER  BOUND  •) 

34  CALL  UNBND 

IF  (L.EQ.1.AND.N0DE0.EQ.2)   GO  TO  46 
IF(L.FQ.l. AND.ABS(R(L) LGT.500)   NODEO  =  2 
GO  TO  8 
41  K=XK 

IF(XMK)  .EQ.2.  )  GO  TO  42 
IF(ABS(XT-X(K ) )-E)  260,210,220 
220  IF(P(L>)  230,8,243 
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240  IF  (XT.LE.X(K))  GO  TO  250 

JS  =  K 

XI (K)  =  +1. 

L  =  L  -  1 

GO  TO  221 
250  ID  =  C 

XLB(K)  =  XT 

GO  TC  42 
210  ID=0 

XU8(K)=XT 

GO  TO  42 
230  IF(XT.GE.X(K) )  GO  TO  210 

JS  =  K 

XI <K)=-1. 

L  =  L  -  1 

GC    TO    221 
260    IF     (R(L) )210, 8,250 
529    WRITE    (6,5701)     L,R(L),Z,ZZ 

WRITE     (  IWR  ITE,5702) 
5702    FORMAT!'     ERROR,    L     IS    NEGATIVE.*) 

46  DO  47  J=1,N 
IF(XSOL( J)  )48, 47,48 

48  WRITE! IWPITE,49) J,XSOL( J) 

49  FORMAT (1HC ,5X,5HXS0L( ,I2,2H)=,F13.4) 

47  CONTINUE 

WRITE  UWRITE,52)ZZ 
WRITE(6,52)ZZ 
52  FCRMAT(1H0,3HZZ=,F1C5) 
528  RETURN 
FND 

SUBROUTINE  AXE 

COMMON  A(  25,  25), B(  25,  25),IB(  25), H<  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB<  25),XI(  25), Y(  25), R(  25) 

COMMON  XSOLl  25),IRV(  25) 

CCMMON  MM,N,M,JS,E,IR,AD,W,Z,IM,DI, I TR , XH , AH , HH, EE , BND 

COMMON  ZZ,DJS, ITE R ,L, I W, NODEO, I NT RES , XT, XK , I  0, RL , BS IGN 

COMMON  IREAD, IWRITE, I8NB 

C  THIS  IS  THE  SUBPROGRAM  WHICH  COMPUTES  THE  A(JS)  COLUMN 
C  NEEDEC  TO  PERFORM  THE  PIVOT  STEP. 

DO  21  1=1, M 
AX(I  )=C. 
DO  21  K=1,MM 

AX (I  )=AX( I  )  +  B( I,K)*A(K,JS) 
21  CONTINUE 

AX(MM+1)=AX(MM  +  1H-A(MM-H,JS) 

RETURN 

END 

SUBPOUTINF  PIVOT 

COMMON  A(  25,  25), B(  25,  25),IB(  25), H(  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB(  25), XK  25), Y(  25), R(  25) 

CCMMON   XSOL(  25),IRV(  25) 

COMMON  MM,N,M,JS,E,IR,AD,W,Z,IM,DI , ITR ,XH, AH , HH, EE . BND 

COMMON  ZZ,DJS,ITER,L,IW,NODEO,INTRES,XT,XK,ID,RL,BSIGN 

CCMMON   IREAD,  IWRITE,  IBNB 

C  THIS  IS  THE  SUBPROGRAM  WHICH  DOES  THE  SIMPLEX  PIVOT  STEP. 

AD  =  AD*ABS<  AX(  IR) ) 
DI  =  AINT(DI*4BS(AX(  IR)  )  +  .5) 
84  WRITE( IWRITE, 214) IR,AX(IR) 
214  FORMAT! 1H0  ,2  5X  ,  3HI R= , I  3, 3X , 3HAX= , F8. 4 ) 
IF(IB( IR) )1,1,2 


36 


K=IB(IR) 

XI(K)=C. 

IB(IR)=JS 

XI (JS)=2. 

E=.5/DI 

DC   2<5    1  =  1,  m 


DC    2Q    1  =  1,  V 

IF (I .EQ.IR)GO    TO    29 

P  =  -AX(  I)/AX( IR) 

H( I)=HU )+P*H(  IR) 

DC    31     J=1,MM 
31    B(  I, J)=B(  I ,J l+P*B(  IR,J) 
29    CCNTINUE 

P=1./AX( IR) 

H(  IR)  =  H(  IR) 

nr     -a  -a       i—  i      u 


SUBROUTINE  VAR 

COMMON  A(  25,  25), B(  25,  25),IB(  25), H(  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB(  25),XI(  25), Y{  25), R(  25) 

COMMON   XSOLI  25),IRV(  25) 

CCMMON  MM,N,M,  JS,E, IR,AD,W,Z, I M, DI , I TR ,XH, AH ,HH, EE ,BND 

COMMON  ZZ,DJS,ITER,L, IW,NODEO, I NT  RES , XT, XK , I D, RL , BS IGN 

COMMON   IREAO, IWRITE,IBNB 

C  THIS  SUBROUTINE  PICKS  OUT  THE  NON-INTEGER  VARIABLE  IN  THE 
C  CURRENT  SOLUTION  THAT  WILL  BE  BOUNDED  AND  BRANCHED  UPON. 

XLDU=0.C 
XLDL=0.0 
IML=0 
IML  =  C 

62  DO    50    IV    =    1,MM 
XDU=1CC0CC0 
XDL=10CCCCC 

IT     ( IB(IV) )    50,50,55 
55    IF(IRV(IB( IV) ) .EQ.0.AND.INTRES.EQ.1)     GO    TO    50 
IF(ABS<X( IB( IV))-AINT(X( IB< IV) )+E) )-E)50,50,60 

60  DO    61    JV  =  1  ,MM 

WRITE    <IWRITE,44)     B( M , JV)  , B( I  V , JV) 
44    FORMAT     (•     B(M, JV)=»,F14.4, «,B( IV,JV)=«,F9.4) 

IF< ABS(B<M, JV) ).LT.E)     GO    TO    61 

IF( ABS(B< IV, JV) ) .LT.E)    GO    TO    61 
IF(B(IV,JV))    64,61,63 

63  XDUJV=ABS(  B(M, JV)/B( IV, JV) ) 
XDLJV=    0.0 

IF(XDUJV.GE.XDU)     GO    TO    61 
XDU=XDUJV 

WRITE     (IWRITE,41)XDLJV,XDUJV, JV,IV,X( IB(IV)),IB(IV) 
GO    TO    61 

64  XDLJV=ABS(B(M, JV)/B( IV, JV) ) 
XDUJV=    0.0 

IF(XDLJV.GE.XDL)    GO    TO    61 
XDL=XDLJV 

WRITE     <IWRITE,41)XDLJV,XDUJV,JV,IV,XUB<  IV))  ,IB(  IV) 

61  CONTINUE 

52    FIDLO    =  X(IBUV))-    A  INT  (  XII  B(  I  V)  )  ) 


37 


53    XLDLIV=xni_*FIDlO 

XLDUIV=XDU*(1.0-FIDLO) 
IF(XLDLIV.LE.XLDL)  GO  TO  66 
XLDL=XLDLIV 
IML=IV 

66  IF(XLDUIV.LE.XLDU)  GO  TO  50 
XL0U=XLDUIV 

IMU=IV 
50  CONTINUE 

IF(XLOU-XLDL)  68,67,67 

67  IM=IMU 
BSIGN=     1.0 
GO    TO    45 

1^0         FORMAT    (Fl.O) 
6  8    IM=IML 

BSIGN=-1.0 
45    WRITE(  IWRITEilO) IM 

WRITE    (6  ,10)IM 

READ    (5,10C)DST0P 

IF    (DSTOP.EQ. 1.0) STOP 

10  FORMATdHC ,  10X  ,     •     FROM    VAR     I  M=     «     ,12) 
WRITE     (6  ,42)     BSIGN 

WRITE  (IWRITE, 40)  F IDLO, XDU, XOL , XLDU, XLOL , IMU, 
1IML, IM 

40  FCRMATd  F  10L0=» , F9. 4 , • , XDU= • F9 .4 , • , X DL= •  , F9 .4, 
,XLDU=#,F9.4, ,,XLOL=« , F9.4, • , IMU=« ,14, 

14,' , IM=» , 14,/) 

41  FORMATC     XDL  J  V=»  ,F9.  4  ,  •  ,  XDUJV=«  ,F9.4  ,  •  ,  JV=  •  ,  14,  •  ,  I  V=«  , 
1I4,,,X(IB(IV))=»,F5.1,,,IR(IV)=«,I4,/) 

42  FORMATC     B  SIGN  =  »  ,F9.4  ) 
IF(  IM.EQ.OJSTOP 
RETURN 

END 

SUBROUTINE  REDO 

CCMMON  A(  25.  25), B(  25,  25),IB(  25), H(  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB(  25),XI(  25), V(  25), R(  25) 

COMMON   XSCL(  25),IRV(  25) 

COMMON  MM,N,M, JS , E , I R , AD, W, Z , I M, DI , ITR ,XH, AH, HH, EE , BND 

CCMMON  ZZ,DJS,ITER,L, IW,NODEO, INTRE S , XT, XK , I D ,RL , BS IGN 

COMMON   IREAD, IWRITE, IBNB 

C  THIS  SUBROUTINE  REMOVES  THE  LAST  BOUND  SO  THAT  6ACK  TRACK- 
C  ING  CAN  START,  AND  RESETS  ALL  OF  THE  OTHER  BOUNDS. 

5  WRITE  (IWRITE, 14) 
14  FORMAT(1HO,10X,  •  FROM  REDO  •) 

DO  21  K=1,K 

XUB(K)=-1. 

XLB(K)=0. 
21  CONTINUE 

LM=L-1 

IF  (L'i.EQ.C)  GO  TO   17 

DO  18  JT=1,LM 

CALL  UNPAK  (JT) 

K=  XK 

IF  (RUT))  19,18,20 

19  XLB(K)=XT 

WRITE(  IWRITE, 11)K,XLB(K) 

11  F0RMAT(1H0,10X,2HK=,I2,5X,4HXLB=,F10.4) 
GO  TO  18 

20  XUB(K)=XT 

WRITE( IWRITE,12)K,XUB(K) 

12  FORMATdH  ,  10X  ,  2HK  = ,  I  2  ,  5X,  4HXUB=  ,  F  10.  4 ) 
18  CONTINUE 

17  JT  =  L 

CALL  UNPAK  (JT) 

ID=500 

If (BSIGN. EQ. -1.1  XT=XT-1 
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IF(BSIGN.EQ.+1.)  XT=XT+1 
CALL  PACK 
BSIGN  =  -  BSIGN 
R(L)  =  -R(L) 
WRITB( IWRITE,15)L,R(L) 
15  FORMAT(lHCf 10Xf2HL=»I2»5X,2HR=fF14.8f • XXXXXXXXXXXXXM 
1  RETURN 
END 

SUBROUTINE  PACK 

COMMON  A(  25,  25), B(  25,  25),IB(  25), H(  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB{  25),XI(  25), Y(  25), R(  25) 

COMMON  XSOL(  25),IRV(  25) 

COMMON  MM,N,M, JS , E , I R, AD , W , Z , IM,OI, ITR,XH, AH, HH, EE , BND 

COMMON  ZZ,DJS,ITER,L,IW,NODEO, I NTRES , XT, XK  ,  I D, RL , BS IGN 

COMMON  IREAD, IWRITE, IBNB 

C  THIS  SUBROUTINE  COMBINES  INTO  A  SINGLE  NUMBER  THE  VARIABLE 
C  BOUNDED,  THE  VALUE  CF  THE  BOUND,  AND  WHETHER  IT  IS  AN 
C  UPPER  OR  LOWER  BOUND. 

IFUT.EQ.O.)  GO  TO  1 

IL=1 

IT  =  XT 
30  IFUT/1C)  11,11,20 
20  IL=IL+1 

IT=IT/10 

GO  TO  30 
11  TL=IL 

RID  =  (  XK*-TL/1C.+XT/10.**(TL  +  1.  )+FLOAT(  ID)  )*(  BSIGN) 

WRITE! IWRITE,10)XT,IL,TL,L,R<L) 
10  FORMATClHOf lCXt  ■  FROM  PACK  ■  , 5X, 3HXT=, Fl C. 4, 5X , 3HIL= 
1,F1C.4,5X,2HL=,I5,5X,2HR=,F14.8) 

GO  TO  2 

1  R(L)  =(  XK  ♦  FLGAT(ID))*(  BSIGN) 

2  RETURN 
ENC 

SUBROUTINE  UNPAK  (JT) 

COMMON  A<  25,  25), B(  25,  25),IB(  25), H(  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB(  25),XI(  25), Y(  25), R(  25) 

COMMON  XSOL(  25),IRV(  25) 

CCMMON  MM,N,M,JS,E,IR,AD,W,Z,IM,DI , I TR, XH, AH.HH, EE ,8ND 

COMMON  ZZ,DJS,ITER,L,IW,NODEO, I NT RES , XT, XK, I D, RL , BS IGN 

COMMON  IREAD, IWRITE, IBNB 

C  THIS  SUBROUTINE  RETREVIES  THE  INFORMATION  COMBINED  IN 
C  SUBROUTINE  PACK. 

RL=ABS(R( JT) ) 
ID=0 

IFIRL.LT.5C0)  GO  TO  10 
RL=RL-500. 
ID=500 
10  IH=RL 
XK=IH 

TL=<RL-XK)*10. 
IL=TL 
RX=IL 

IXT=(TL-RX )*10.**RX+.5 
XT=IXT 

BSIGN  =  SIGN(1.0,R(L>  ) 
RETURN 
END 
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SUBROUTINE  BOUND 

COMMON  A(  25.  25). B(  25,  25) • I B(  25), H<  25),AX(  25) 

COMMON  X(  25),XLBt  25),XUB(  25), XK  25), Y(  25), R(  25) 

COMMON   XSOL(  25),IRV(  25) 

COMMON  MM,N,M,JS,E,IR,AD,W,Z,IM,DI , I TR , XH , AH , HH , FE , BND 

COMMON  ZZ,DJS,ITER,L,  IW,NOOEO, INTRES , XT , XK  ,  I D ,RL, BSIGN 

COMMON   IREAD, IWRITE,IBNB 

C  THIS  SUBROUTINE  SETS  THE  VALUE  OF  BOUND  OF  THE  NEXT  BRANCH 

IDBND  =  ID 

LSET  =  C 

IF(  ID.EQ.500)  GO  TO  2 

L  =  L+1 

LSET  =  1 
23  CALL  VAR 

IF  (  IM.GT.O)  GO  TO  17 

WRITE  <IWRITE,14)  IM 

L  =  L  -  1 

RETURN 
17  K=IB(IM) 

Ih=X(K)*Dl+.5 

XT=IH 

X(K)=XT/DI 

IH=X(K)+E 
32  XT=IH+1 

IF(BSIGN.EQ.-L)  GO  TO  31 

CALL  PACK 

GO  TO  6 
31  Y( IM)=XT-X(K) 

X(K)=XT 

XLB(K)=XT 

IB(IM)=-IM 

M=M+1 

XI(K)=0. 

DO  1  J=1.MM 

B(M, J)=Bl IM, J) 

1  B( IM,  J)  =  -BUM,  J) 
H( IM)=-H( IM) 

W  =  Y(  IM) 

IW=0 

WRITEUWRITE,16)  W 

IF ( IDBND. EQ. 500)  RETURN 

XK=K 

CALL  PACK 

GO  TO  3 

2  K=  XK 

IF(ABS(XT-X(K) )-E)  21,21,22 

21  ID=C 
IF(BSIGN.LT.O.O)  GO  TO  30 
XUB(K)=XT 

GO  TO  23 
30  XLB(K)=XT 
GO  TO  23 

22  IF(XT.GE.X(K) .AND. BSIGN. EQ.+l.)  GO  TO  21 
IF(XT.LE.X(K) .AND. BSIGN. EQ.-l. )  GO  TO  21 
DO  5  1=1, M 

IF(IB( I) .EC.K)GO  TO  4 
GO  TO  5 

4  IM  =  I 

WRITEt I  WRITE, 14) IM 

14  FORMAT(1HO,10X,    ■     FROM    BOUND    I M=    ■     ,15) 
GO    TO    8 

5  CONTINUE 
WRITE( IWRITE.15) 

15  FORMATdHClOX,     ■     ERROR    •     ) 
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8  ID  =  0 

IF  (BSIGN.LT.O.O)  GO  TO  31 

6  Y( IM)=XIK)-XT 
X(K)=XT 
XUB(K)=XT 
WRITEIIWRITE,10)K,XIK) 

10  FORMAT(1HO,10X,  ■  FROM  BOUND  ■  , 5X ,2HK= , I  2  ,5X  ,2HX= , F10 
IB(IM)=-IM 
M=M-U 
XI  (K)  =  0. 
DO  7  J=1,MM 

7  BIM, J)=-B( IM,J ) 
W=Y( IM) 

IW=0 

WRITE(  IWRITE, 16)W 
16  FORMAT!  1HC10X,  •  FROM  BOUND  W=  •  ,F10.8) 
IFILSET.F0.1 )RETURN 
L  -  L  ♦  1 
XK  =  K 
CALL  PACK 
3  RETURN 
END 

SUBROUTINE  UNBND 

COMMON  A(  25.  25), B!  25,  25), IB!  25), HI  25),AXl  25) 

COMMON  XI  25),XLB(  25),XUBI  25), XK  25), Y(  25), R{  25) 

COMMON  XSOLI  25), IRVI  25) 

COMMON  MM,N,M, JS , E , I R , AD , W , Z , I M , D I , I TR.XH, AH ,HH, EE , BND 

COMMON  ZZ,DJS,ITER,L, IW,NODEO, INTRES , XT , XK , I D,RL , BSIGN 

COMMON  IPEAD, IWRITE, IBNB 

C  THIS  SUBROUTINE  RESETS  ALL  THE  BOUNDS  AFTER  EACH  BACKTRACK 
C  AND  DETERMINES  WHICH  NODE  TO  BACKTRACK  TO. 

6  IF(ABS(R(L))-500.0)  1,1,2 

1  CALL  REDO 
ID=5C0 
RETURN 

2  IFIL.EQ.l)  CALL  REDO 
IF(L.EG.l)   RETURN 

L  =  L  -  1 
GO  TO  6 
ENC 

SUBROUTINE  ELEM 

COMMON  A(  25,  25), B(  25,  25),IB(  25), HI  25),AXI  25) 

COMMON  XI  25),XLB(  25),XUB(  25), XK  25), Y(  25), R(  25) 

COMMON  XSOLI  25),  IRVI  25) 

COMMON  MM,N,M, JS,E,  IR,AD,W,Z,  I M, DI , ITR ,XH, AH, HH, EE, BND 

COMMON  ZZ,DJS,ITER,L,IW,NODEO,INTRES,XT,XK,lD,RL,BSIGN 

COMMON  IREAD, IWRITE, IBNB 

C  THIS  IS  THE  SUBPROGRAM  WHICH  OBTAINS  THE  PIVOT  ELEMENT 

C  FOR  THE  GIVEN  COLUMN  DETERMINED  EITHER  IN  'LP*  OR 

C  IN  'BOUND*. 

C  IT  GIVES  WHICH  VARIABLE  WILL  LEAVE  THE  BASIS. 

IR=C 

CALL    AXE 
EE=E 

WRITE! IWRITE, 43)E 
43    FORMAT(1HO,10X,12HFROM    ELEM    E=,F10.4) 
IFIBND.EQ.C. )G0    TO    12 
IFIXK  JS).EQ.1.)G0    TO    1 
GO    TO    2 
1    XH=XUB!JS)-X(JS) 
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IF(XUBUS)  .EQ.-l.  )XH=  1000000. 
GO  TO  12 

2  XH=X( JS)-XLB( JS) 
12  DO  25  1=1, MM 

IFIBNC.GT.C. )G0  TO  810 

IF(AX( I) .LE.E)GO  TO  25 

AH=AX(  I) 

HH=H(I  ) 

IF(IR)25, 811,812 
311  IR=I 

XH=HH/AH 

EE=E/AH 

GO  TO  25 
812  CALL  MIN  (I) 

GO  TO  25 
810  IF(ABS(AX( I) ).LT.E)GO  TO  25 

IF (X  I ( JSJ.EQ.I.IGO  TO  8 
C  DECREASE. 

IFUB(I)  )3,25,4 

4  K=IB(I ) 

IF(AX( I) )5,25,6 
6  IF(XUB(K). EQ.-l. ) GO  TO  25 
HH=XUB(K)-X(K) 
AH=-AX( l)*XI< JS) 
CALL  MIN  (I) 
GO  TO  25 

5  HH=X(K)-XLB(K) 
AH=AX(  I)*Xl( JS) 
CALL  MIN  (I) 

GO  TO  25 

3  IF(AX( I) 115,25,16 
16  HH  =  Y(I  ) 

AH=-AX( I)*XI( JS) 

CALL  MIN  (  I) 

GO  TO  25 
15  HH=Y(I  ) 

AH=AX(  I)*XI( JS) 

CALL  MIN  (I) 

GC  TO  25 
C  INCREASE. 

8  IFUBI  I)  )30,25,4C 
40  K=IB(I ) 

IF(AX( I) )6,25,5 
30  IF(AX(  I)  )16, 25,15 
25  CONTINUE 

RETURN 

END 

SUBROUTINE  MIN  (  I  ) 

CCMMON  A(  25.  25). B(  25,  25),IB(  25), H(  25),AX(  25) 

COMMON  X(  25),XLB(  25),XUB(  25),XI{  25), Y(  25), R(  25) 

CCMMON  XSCL(  25),IRV(  25) 

COMMON  MM,N,M, JS, E , IR, AD , W, Z , I M, 01 , I TR ,XH, AH ,HH, EE , BND 

COMMON  ZZ,CJS, ITER,L, IW,NODEO,  INTRES, XT, XK , I D,RL , BSIGN 

CCMMON  IREAD,IWRITE, IBNB 

C  THIS  IS  THE  SUBROUTINE  USED  TO  CALCULATE  THE  MINIMUM  THETA 

IF (ABS(HH-AH*XH)-EE) 70,70,2 

2  IF(HH-AH*XH)24,70,25 
24  XH=HH/AH 

3  IR  =  I 
EE=E/AH 
GC  TO  25 

7C  IFUR.EQ.O)  GO  TO  25 
IF(IB( I))74,25,73 

73  IFUR.EQ.O )  GO  TO  25 
IF  (  IB (IP)  )  25,2  5,4 

74  IFIIR.EQ.O)  GO  TO  3 
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IF    (  IBUR)  )    4,25,3 
4    IF{ABS(AX(IR)  >-ABS(AX(I))  )  25,  2  5,  3 

25    RFTURN 
END 
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ASSEMBLY  LANGUAGE  TIDING  ROUTINE 


SETIME 


FIXUP 
fcATIME 


SAVE 
TEMP 
TIME 


( 14, 12) 
12,15 
SETIME,12 
13, TEMP 
13, SAVE 

TASK,FIXUP,TLINTVL=TIME 
13, TEMP 

T 


START 

SAVE 

LR 

USING 

ST 

LA 

STIMER 

L 

RETURN  (14,12) 

DROP   12 

SAVE   (14,12) 

RETURN  (14,12) 

ENTRY  WATIME 

SAVE   (14,12) 

LR     12,15 

USING  WATIME,12 

L      2,0(1) 

ST     13, TEMP 

LA     13, SAVE 

TTIMER  CANCEL 

L      1,TIME 

SR     1,0 

ST     1,0(2) 

L      13, TEMP 

RETURN  (14,12),T 

DS      18F 

DS      F 

dc    x'osocoooc 

END 
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